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METHOD AND APPARATUS FOR MEASURING LOCALLY AND SUPERfl- 
CIALLY THE SCATTERING AND ABSORPTION PROPERTIES OF TURBID 
MEDIA 



)h«5 option isajTUf icr/dMijoolHf 

1 . Field of the invention 

The present invention relates to a method and an apparatus to quantify the optical 
scattering and absorption properties of turbid media, that can be applied in the extended 
optical domain of electromagnetic waves from far infrared (up to the microwave domain) 
to the far UV domain, i.e. where the energies of the photons composing the 
electromagnetic waves range from a few milli-electron-Volt to more than ten electron- 
Volt. More precisely the present invention relates to a non-invasive measurement, over a 
small area of the sample surface, Local and superficial characterization of biological 
tissues Ln vivo is a major application of this invention. 



2. Related Background Art. 



Different techniques have already been proposed to quantitatively determine the 
absorption and reduced scattering coefficients of turbid media 1 . Most of the non- invasive 
methods are based ^on the measurement of spatially and/or temporally-resolved 
reflectance. The principle^SNas follows: the turbid medium is illuminated by a collimated 

ft. ' — — '■' ' " I ... " " ~ - , 



or focused li ght source .v ui &backs ra ftered light is measured fay one or several detectors 

— -v^ " : 



Different types of measurements are possible, depending on the time-dependence of the 

\ 

illuminating source: steady-state (continuous source), time-domain (short pulsed source) 
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or frequency domain (amplitude modulated source)^The.prej»eftt invention relates to the 
case of steady-state measurements, performed at^dififerenc distances between the source 
and the detectors. However, the technique presented here can be complemented by tinse- 
or frequency-doma^measurements. n ^ ^ / *f- _ ^ f - ^ ^ _Sf^ 




The range^of p values is an important point to consider, when comparing different 
methods based on the measurement of the reflectance. First, the probed volume of the 

hi.hu gwii.il Hi imintmruM, 

turbid medium is related to tjie source-detector separation p. (The larger the source- 



detector separation, the deeper the average depth probed Second, depending on the range 
of p, different mathematical processing must be used to obtain the optical properties from 
the raw data. 
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At least two cases must be distinguished. 



I) The first case corresponds to source-detector separations larger than several transport 
mean free paths. Fo^ typical biological tissue optical properties 2 , this case correspond to a 
source-detector separation larger than 2 mm. An analytical form of the reflectance can be 
obtained^ro^^ 

than th^ scattenng coefficient p., ^typically ten times). In such a case, the relevant optical 



propertiesjfe the^^aSiveln^ the absorption coefficient and the reduced scattering ^ JL 
:oefficient. The average depth of probing is on the same order than the source-detector ' 



separation p. 
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Such methods have been already published, and are the object of patents^Ref , Patent 




5 T 5 17,987 Tsuchiya, Patent 5,676,142 Miwa et al.). 



2) The second case corresponds to source-detector separations close to one transport 



mean free path. For biological tissues*, such source-detector separations correspond 
typically to distances from 0.1 to 2 mm. The average depth probed is on the order of 1 
nun Such small source-detector separations enable the measurement of the optical 
properties of a small tissue volume. 

Wang et al. (Patent 5,630,423) proposed a method for the determination of the reduced 



scattering coefficient only, using an optical beam of oblique incidence. Moreover, their 
analysis does not include the effect of the phase function. K essler et al. (Pat gotS- 
5,284,137 and 5,645,06 Ij jsroposed a method for measuring the local dye concentration 
and scattering parameters in animal and human tissues, based on spatial and spectral 
measurements. However, their methods do not enable the simultaneous determination of 
the absorption coefficient, reduced scattering coefficient. Other publications may be of 



1 . See "Welch, A. J.; van Gemert, M. J. C. Optical Thermal Response of Laser Irradiated 
Tissue; Plenum publishing Corp., New York, 1995", and references therein. 

2. W.-F. Cheong, S. A. Prahl, and A J. Welch, "A Review of the Optical Properties of 
Biological Tissues," IEEE J. Quantum Electron. 26,2166-2185 (1990). 




concern; 
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3. Summary of the Invention: /V\ & n ( > I* 

In the case of source-detecror separations close to one transport mean free path,_light 
propagation has been mod el ed using Monte Carlo simulations. In contrast to the previous 
case where the source-detector separation is larger than several iransport mean free paths, 
it was discovered by the inventors that not only the first moment of the phase function 



muse be considered, but also the second moment More precisely, the relevant, optical. 





— JC properties are the refractive indexi^^bsorptio^^^r^uced scatterin g coefficient^an d a 
(^new parameter y^l-g2>/(l^g^ takes into ac count the first Wo m oments,of.the^phase 



function, denoted gi and %i respectively. 
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Different methods have been proposed for local characterization of turbid media, and in 
particular of biological tissues. Nevertheless none allows for the simultaneous 
determination of the absorption, the reduced scattering coefficient and the parameter y. 



In the present invention we present a method and apparatus for the measurement of the 
absorption coefficient, the reduced scattering coefficient and the said phase function 



parameter 7, from the spatially-resolved : reflectance data at short source detector 



separation. These three parameters, which can be measured at different wavelengths, 
enable us to characterize turbid media, such as biological tissues. Our method is based on 
a theoretical study of the reflectance at short source-detector separations, that we 
performed vji^Mcote Carlo simulations. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 , Description of the principle of spatially-resolved reflectance measurement. 

p denotes the source - detector distance, CD fi0 wce and eodeteoor are ^ e solid angles of light 
acceptance determined by the Numerical Apertures (N.A.) of the optical fibers 



Fig.2 Probability density function of the mean depth of scattering events per detected 
photons. 
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i 

IB 

y Fig. 3. Examples of spatially resolved reflectance curves obtained with different phase 



functions. Case of matched refractive index (n=1.0). 



Fig. 4. Monte Carlo simulations and fits of the spatially resolved reflectance curves by 
mathematical expression in the formi?(/?) ~^A^p i /r Jt y^// f B(/4 a )J . Case of constant 

;=1.9 and 1 mm -1 , numerical aperture of the source and 

detectors = 0 37. A - 0.0647/^ J24 exp(-0J61/>), B * 0.18653 0.8466 pj + 1.836/t,' 

Fig.5.a. Basic description of the apparatus: Case of measurements with an optical fiber 
probe. The apparatus comprises a probe connected to a light source unit and a detection 
unit, both units being connected to a computer or signal processing unit The optical fiber 
probe is constituted of one or more illuminating fibers, and one oi more collecting fibers. 
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Fig.5.b.Ba$ic description of the apparatus: Particular embodiment where the probe is 
composed of n delecting fibers (typically 6) in parallel. A broadband or white light source 
is used for excitation and the retro-diffused light is collected by the n detecting fibers in 
parallel and dispersed in parallel in a spectrograph, before being detected on a 2D CCD 
camera and transmitted numerically to a PC for xm/> m* determination at each 
wavelength Z . 



*fi' Fig.S.c. Basic descriptiou of the apparatus. Case of a measurement with sources and 

m 

Q detectors directly in contact with the turbid medium. 

m 

- = * 

li 

□ Fig.5.d. Basic description of the apparatus. Case of non-contact measurements with a ID 

Q or 2D detector, coupled with an imaging device. Scanning systems 1 and 2 are optional 

m 

O they may be distinct from one another or confounded in a single scanner. 



Fig.5,e. Basic description of the apparatus. Case of non-contact measurements with 
source and detection beams deflected by one or two scanning devices which may be 
distinct from one another or confounded. 



Fig.6. Examples of the sample side of the fiber optical probe, 

6a. Simple arrangement for a single measurement 
6b. Arrangement for symmetrical measurements. 
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6c. Arrangement for multiple measurements. 

6d. Arrangement for multiple measurements, using an multicore fiber, or an optical fiber 
bundle. 

6c. Arrangement for symmetrical measurements where detecting fibers are arranged on 
an elliptic path so as giving a regular spacing from the illuminating fibers to the detecting 
fibers, while maximizing the collecting surface of the fibers. Additional fibers F can be 
put in the center of the probe or around the probe for other light channels: collecting 
fibers for fluorescence or Raman scattering for example. 

Fig.7. Example of a measurement of the spatially-resolved reflectance with the probe 
shown in Fig.6.a. The measurement, performed on a microsphere suspension, is 
superimposed to a Monte Carlo simulation. The optical properties, computed from 
published water optical properties Cu n ) and Mie theory (ju s ', are: n=L33, /< ff =0.0004 

m 

G mm* 1 , ( u/=1.0 mm" 1 , ^=2.2. The calibration was performed with a siloxane sample of 

known optical properties. 
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Fig.8. Relation between the parameters R(p-1 mm) and \d pinR(p~l mm)\ and the 
optical coefficients and fx Q .Case of ^==1.5 and 1.9. Probe of refractive index = 1.5, 
sample of refractive index - 1.4, optical fibers diameter = 200 jim, ^=0.37 (source and 
collection). 
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Fig.9. Plot of the parameter ~->/R for different /values (L0 t 1.8, 2.5), different reduced 

dp 

albedo a* (1, 0.95, 0.9) and for fixed /?=1 ram. Mismatched refractive index n=*IA 



Fig. 10. Plot of yjR (a ' = 1} -^(a') for different /values and reduced albedo a\ and for 
fixed p=l mm. Mismatched refractive index n=1.4. 

DETAILED DESCRIPTION OF THE INVENTION 



1. Definitions 

The presenttdefinitions are given for clarity of the concepts developed in the patent and 
are wordecHE agreement with ref. 1 . 



The concept of spatially resolved reflectance is illustrated in Fig.l. Consider a. collirnated 
light beam impinging on a turbid medium (through air or through a light guide, or even 
generated locally), in a given solid angle <b 50utoc determined by the Numerical Aperture of 
the illuminating fiber: N.A. which is defined as the half angle of the entrance angle in the 
fiber, and on a limited area A sour{ x> where the source irradiate the surface. The spatially 
resolved reflectance R(p) is defined by the backscattered light power in a given solid 
angle a>deiecier* at a distance p from the source, per unit area and normalized by the source 
power. The collecting area of the detector or detecting fiber is A<t 4 u* tCr . The distance p is 
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referred as the source-detector separation (also in the case when light guides or imaging 
devices are used for the illumination and for the collection of the backscattered light). 

R(fi) depends on the optical properties of the turbid medium, defined below Note that 
R{p) depends also on the source and detectors characteristics, i.e, numerical aperture and 
sizes. In a general situation, R($ is the result of the convoluri on_^fLfcjl_j^jTe^ 
.reflgctMce ^urve and the source irradia nce and detector sensitivity over, respectively, the 
surfaces A SOyrcc and A^ tlttlor . When A soltrce and Adttectw are smaller than the required spatial 
resolution , R(p) is confounded with the corrected spatially resolved reflectance. 



It is commonly admitted that the fundamental optical properties of a turbidnie 



determined by th^av erage index of refr action n of the mediun^ ^^^^^^I^^i^L 
^^^^^^^^^^^^^iid tl^^^^fti ^tion j7f ^where & is the scattering 
angle. The absorption coefficient p a [mm" 1 ] is defined as the probability of absorption per 



unit infinitesimal pathlength. The scattering coefficient fi s [mm* 1 ] is the scattering proba- 
bility per -unit infinitesimal pathlength. The phase function p($ is the density probability 
function for the scattering angle 0. The phasefi inction is nor malized as follows: 



p(0)s\nffd# 



(1.1) 



In turbid media and in biological tissues in particular, the phase function p(&j\z a 
continuous function of the scattering angle 8 and generally depends on the dielectric 



properties and on the material microstructure. It can be approximated by a development 
in Legendre polynomials: 
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Where jg„ is the order mome nt of t he phas e function and P n (Q is the Legendre 
polynomial of order n: 



(13) 



HI 



The first moment of the phase function is also called thefanisotropy factor, and is often 
simply noted g ( & gj). It represents the mean cosine of the scattering an'gH The reduced 




scattering coefficient ju s ' is aenned asr * " 

The transport mean fires path (or reduced mean free pat h) mfjp * of the photons, in a tuibid 



medium . is defined as: 



(1.5) 




The reduced albedo a ' is the ratio: 



(1.6) 



As a result of the present invention described in more details in the next section, it is also 



necessary to define anothdrphase function para meter c alled n 





(1.7) 
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All che parameters listed above are referred as optical properties. They are wavelength 
dependent* and can vary in space and time! V l2TEe^foU^^ gj and g* will be 
considered as parameters characterizing the turbid medium. 




Si 



m 



2. Reflectance measurements at distances close to one transport mean free path 



V 



The methods described in the US Pate^T^ are based on 

measurements wim large range of source-detector separations, typically from I mfp f to 
10 mfp\ tn such cases, the volume probed is on the order of 10-1000 (mfp*) 3 . In contrast 



with such a large so^le investigation, the present invention is directed to a novel approach 
where the volume probed is much smaller, on the order of 1 (mfp 1 ) 3 . This is achieved by 
using only smailj^frce detector separations, typically from 0.1 to 2 mfp\ The lateral 
dimension of probini is limited to this range of distances. 



Note that the present invention can be complemented by large source-detector separation 
measurements, if the optical properties of both superficial and deep parts of a turbid 
medium must be determined. 



A model of photon migration in tissue was necessary to predict the relationship between 
the measured reflectance and the optical properties. Analytical solutions from the 
diffusion equation are not appropriate in our case because we are interested in the 
reflectance close to the sourc^a^distances comparable to the transport mean free path 
[mfp*]. We have performed Monte Carlo simulations to predict the measured reflectance 
of an homogeneous semt-infinite-tufbidrmedia". 
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The exact diameter of the illuminating and collecting fibers, as well as their Numerical 
Apertures (N.A.)> have been taken into account in the simulations. The mismatch of 
index of refraction at the surface of the medium have been also taken into account in the 
simulations, by using the Fresnel law for each photon reaching the surface. 



This is also one result of this simulation to compute the average depth of probing, 
illustrated in Fig.2 A It is demonstrated that only the superficial part of the turbid medium 
is probed if small source-detector separations are used For this, we determined the depth 
below the surface ofleach scattering event in the simulation. With this information we 
determined the average depth of all the scattering events for each detected photon, which 



we present as a probability density function in Fig.2. . This figure shows that the average 
depth of scattenng^^proximately around 1 mfp\ Moreover it showed that the part 



located below 2 mfp* were not playing a significant role in the measured signal (for p< 
1.5 ixifp'). 



The spatially resolved reflectance RQjij, w ith short source-detector separations is' more 
complex than in the case of large source-detector separations. Indeed, for small source 
detector-separation conditions, the inverse problem, i.e. calculating the localized 
absorption and reduced scattering coefficients, is sensitive to the scattering phase 
function. It is part of the present invention to have shown, from Monte Carlo simulations, 
that only the first and s econd m oments of the phase function can be taken into account. 
Moreover, it was established that the influence of these two moments are not 
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independent Indeed, it is a merit of this invention to demonstrate that only one parameter 
?^fl-g2)/n-gi}tyhich depends on the first and second moment* of the phase function 



can be used to characterize accurately the reflectance at short source-detector separation. 




FigT3"iliuitrates tne fact thaTthe parameter / is the only predominant parameter of the 
phase function thattaust be taken into account (and not g=*gj as frequently mentioned in 
literature). Different* reflectance curves, obtained from Monte Carlo simulations are 
shown in Fig.3. Four aifferent phase functions were used for the simulations. Three phase 
functions are characterized by ^1-25, but different gi and g 2 values, Fig.3. shows that 



almost identical reflectance curve are obtained for distances pfi x 1 >0.3, if the gj and £?are 
varied in such a way that* the parameter ;k stays constant (^=1.25). For comparison, the 
reflectance computed ^&h$a different g(~g<) value is also presented (r=2.25). It appears 
then clearly that the parameter g ( =grf plays a significant role in the reflectance. 



-Tie-pararaeter— ;<4epends-on~&^ of the set of scatterers inside the turbid 



medium: shape, distribution of sizes and refractive index. 



The possibility of determining rhis parameter ^is an important part of this invention, and 
has never been reported earlier. The deter mination of r can be related-to-the mieiywgripi^ 
structure of the sample, since y is related to the sizes and refractive indexes of the 
scatterers. 
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In most cases, die average refractive indices n and die source and detectors characteristics 
are fixed and known. In such case, only three parameters are determined from measure* 
ments at short source-detector separation: f*s ' and /. 

It is important to note that the parameter y y cannot be estimated if large source detector 
separations 3re used. Indeed, the reflectance is only sensitive to ^at short source-detector 
separation. 

a 
in 

-ijgs It waS'alseiierived-fromJ^gnte Carlo simulations that the reflectance/?^, for/?* lmfp', 

m 
m 
w 



D 

Q 

m 

Q 



can be reasonably well approximated by the expression: 

R(/>)s[A(p,K,r)+B(#..,n)Y (1.8) 

In many cases, we found that Equ.( 1.8) can also be written as: 

R(p)s[A(p,K,r)+HB(M.)f (1-9) 

it is important for the description of the invention to evidence the fact that the function A 
depends on & the scattering properties (i.e. and j) but not on /x fl . In contrast, the func- 
tion B depends on t u a and u s ' but neither on rnor p. An example of Equ. (1.9) is shown in 
Fig. 4. 

This formulation is of great help to solve the inverse problem, as described in section 5.3. 
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3. Apparatus 



In the first embodiment, The apparatus can be divided in three parts, described in Fig.S.a. 



in 

ifl wavelengths. 



The first part islthe illuminating system. Any light source can be used. For examples: a) 

white sources, such as halogen or xenon lights, metal halides or fluorescent or 
\ 

phosphorescent sources, b) sources such as lasers, laser diodes, optical fiber lasers, light 
emitting diodes^p^perluminescent diodes, c) the sources described in point a) and b) 
where monochromators, filters or interference filters are added to select a given set of 




3 In the first preferred embodiment, The light power is conducted to the investigated 



sample by the probe, w^ich is the second part of the apparatus. The probe is preferably 



made of optical fibers, tc^ illuminate and to collect the backscattered light. But GRIN rods 
or other type of light pip^s can also be used. Different possible arrangements of optical 
fibers are illustrated in Fig^6 . Two different modes of measurements can be chosen. First, 
one fiber is used to iHumiifate the sample and at least two other fibers are used to collect 




the backscattered Hghtjft mp different distances. Second, one fiber is used to collect the 
light and at least two other fibers, located at two different distances from the first one, are 
used to illuminate sequentially the sample. 



The arrangement of the different fibers can be replaced by any imaging system or image 
guide, such as multicore optical fibers, or optical fiber bundle. 
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The light collectedly the probe is analyzed by the detection unit, which is the third part 
of the apparatus. rf wide spectral light sources are used (such as halogen or xenon lights), 
a spectrograph can fee put between the probe and the detector to get wavelength 
dependence of the badcscattered signal (either Ln the source or detection unit). Any types 
Q. of detectors can be Used. For example, photodiodes, avalanche photodiodes or 

photomultiphers can b&jpsigned to each collecting fibers. Simultaneous detection of each 
12 collecting fiber can aisojbe achieved using linear or two-dimensional detectors such as 

Charge-Coupled Detectojs (ID or 2D), intensified CCD or array of photodiodes. 



A particular embodiment is described in Fig.5.b. : the probe is composed of n detecting 
fibers (typically 6). A broadband or while light source is used for excitation and the retro- 
diffused light is collected by the n detecting fibers and dispersed in parallel in a 
spectrograph, before being detected by a 2D CCD camera and transmitted numerically to 
a PC for y,ji z fx a determination at each wavelength k. This particular embodiment has 
the advantage of yielding the wavelength dependence of y, H*', fi a from a single 
measurement. 



A second type of embodiment is presented on Fig. 5.c. The difference with embodiments 
presented in Fig. 5.a. and Fig. 5.b. is that optionally no optical fibers, light pipes or grin 
rods is used. The light source unit is directly in contact with the turbid medium, as well as 
the detector unit. Com&fa%f optics, micro-optics or imaging optics (DOE Diffractive 
Optical Elements for example can be put between the turbid medium and the actual light 
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sources and detectors. The different type of sources and detectors cited in example for the 
first embodiment can be i^^chfor the second type of embodiments. Hybrid design, such as 
arrangements mvolvih/b^mdirect contacts sensors or detectors and fibers, light pipes or 



grin rods are also included 



in the present embodiment. 



The third embodiment is described in Fig.S.A and Fig.5.e. The only difference with the 
first and second embodiment is that non-contact measurements are performed A 
coilirnating system allows for a point-like illumination on the turbid medium. An imaging 
system enables the measurement of the spatial distribution of the reflectance. The 
detectors can be either aiWray (ID or 2D) of detectors (Fig.S-d) with optional scanning 
systems 1 and 2 which cany operate separately or can be confounded in a single scanner, 



or a single detector (Fig.S.e). In this last case, a scanning device is used to obtain the 
spatially -resolved reflectancel|A fiber bundle, multicore fiber or relay optics (grin rod or 
multiple lenses) can be pUKs|^een the focal point of the imaging system and the 
detectors). The different^qggJof sources and detectors cited in example for the first 
embodiment can be used for the third embodiment 



X 



Spatial images \>f the parameters pi/ and y can be obtained by a series of 
measurements atXdifferent locations, that we call multi-site measurements. Some 
mechanical or opticdpx^nning device can be used for this purpose. The resolution of 




such images is oh theWder of the mean source-detector separation used for a single site 
measurement. AlliBSjje embodiments can be expanded to perform multi-site 
measurements, by duplicating and/or multiplexing the illuminating or measuring devices. 



\ 
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For example, the optical fflbgr probes shown in Fig.6., can be duplicated and put side by 
side. The scanning systemYescribed in the third embodiment can also be expanded to 
perform multi-site measurements. 

4. Normalization and calibration 

The differences of transmission between each fiber are corrected by performing a mea- 
surement on a turbid phantom illuminated uniformly. 

The background light, measured with the light source turned off, must be subtracted from 
the signal. 

In order to perform absolute intensity measurements, a calibration is performed on a 
turbid medium of known optical properties. Examples of such a medium are: 1) solid or 
liquid turbid medium which properties have been measured by other techniques, or 
reported in the literature, 2) water suspension of microsphere of known size distribution 
and refractive index. Absorbing dye may be added to the suspension. In case 2) the 
scattering properties are calculated from Mie theory, and the absorption coefficient is 
assumed to be equal to the water absorption coefficient, if no absorbing dye is added. If 
an absorbing dye is used, the absorption coefficient can be measured by a 
spectrophotometer, before mixing the solution with any scattering materials. 

A Monte Carlo simulation is performed with the optical properties of the calibration 
sample. The simulation is then divided by the experimental reflectance performed on the 

18 



AMENDED SHEET 



7- 1 1 -2000 CH 009900476 




O 



calibration sample. The result, that must be independent of the source-detector separation, 
is defined as the calibration factor. Each new measurement is multiplied by the 
calibration factor. 



An example of a measurement, after calibration, is shown in Fig, 7. The measurement was 
performed with the apparatus described in the first embodiment, with a probe similar to 
the one described h Fig.6.a. A prior calibration was performed on a solid turbid medium, 
which optical properties were measured by another technique (frequency domain 
measurement, cf, ReA .). The sample was a water suspension of polystyrene microsphere 



of 1 \rrn diameter. The measurement is compared to a Monte Carlo simulation performed 
Y with the scattering properties calculated from Mte theory and absorption coefficient of 



water, shows an eX$<emqitt agreement between the measurement and the simulation, 



IJl calibration. 



which confirm the validitjjj of our simulation model, experimental measurement and 



5. Signal processing 
5.1. Contrpl of the homogeneity of the area probed 



Artifacts during a ^measurement, for example due to bad contact between the probe and 
the sample, or heterogeneity of the sample, can be detected by the following optional 
procedure. Two illuminating fibers are disposed s)Tnmetrically in regard to the collecting 
fibers (see Fig.6.b). I^/the\ sample is homogeneous, the reflectance curve should be 
y identical with eitner/inljnimating fiber. Therefore, heterogeneity of the investigated 

region or obstructions beneath the fibers are detected by comparing the two curves. If the 
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<&y^ two curves arc suffirieAjjjy^ose, the measurement is validated and the average of the two 
^jv3i^ V curves is calculated. 



m 

Id 



m 

.'KT. 

M 



5.2. Smoothing procedure and computation of the derivative of Ac curve . 

Functions in the form m x p m% exp(mj/?) were always found to fit well Monte Carlo 
simulations for a restricted range of distances. Therefore, a smoothing of the 
experimental reflectance R(p) can be processed by fitting R(p) to m^p"* exp(/w 3 /?) . 

The determination of the slope of the logarithm — (jxxR(p,/f f7 a is also derived 
from this fit, using the following formula; 

— (bA(A/fJ.fl-,r))-^- + m, (1.10) 



The slope of the square root of R(p) can also be derived (the use of these slopes are 
described below in section 5.3 ): 



_d_ 

dp 



(1.11) 



5.3. Inverse problem 

The so-called "direct problem" consists in computing the spatially resolved reflectance 
R(p) from the values of the optical coefficients n. ft a ,, fi s '. y involved in a model of 
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propagatioa of the light in turbid medium and whereby the said "model" incorporates a 
Legendre polynomial development to the second order of the said "phase function", and 
whereby the said "phase parameter" y is introduced in the computation as an 
independent parameter. The so-called "inverse problem" consists in extracting the optical 
coefficients n, /i Cft Mj\ y from the spatially resolved reflectance data R(f$ . 

We developed different methods to solve the inverse problem, which consists in 
extracting optical coefficients from the reflectance data. 

Method 1. 

Monte Carlo simulations show that the measurements of the reflectance intensity R(p) 
and the slope of \nR(p) (denoted dplnR\ determined at a fixed distance p, can be used 
to derive and Ma for a given /value. Fig.8. shows graphically the relationship between 
M/ and ju a and the two parameters R(p=l mm) and \dplnR(p=t mm)\ in the case of 
p=l mm, y- 1.5 and 1.9. Note that any other choice of p is possible, as long as p is on 
the order of I mfp r . We see in Fig. 8. that ^ 9 and /i a can not be determined uniquely if y 
is unknown. Further insight for optimizing the inversion strategy is provided by three 
additional features of Fig.8. 



First, the determination of fi s ' is only weakly influenced by y Indeed in Fig.8. the errors 
induced by error in /are typically 10% for /*/ around 1 mm' 1 . Second, although absolute 
determination of \i h is not possible when y is not precisely known, relative absorption 
changes can be still evaluated. Third, the mdetermination of the parameter y may be 
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resolved by the values of R(p) and/or j dplnR(p)\ at other distances. Therefore the follow- 
ing procedure can be used: 

(1) determination of/;,' and^ a from R(p~l mm) and \d plnR(p~l mm)\ forasetof 
values / For example: r=l.O, 1.1, 1.2,... ,2.5. 

(2) simulations with the different sets of fi s ' and p a obtained from point 1) 

(3) comparison between the simulations and the reflectance profile for distances 
0.3</><2 mm. 

This iast step allows us to determine the correct values of y, p s ' and u*. Points 1 to 3 can 
be done iteratively to evaluate rmore precisely, using a finer discrimination of /values. 

This method described by the iteration of points 1) to 3) corresponds to finding a 
simulation that fit closely a measured curve R(p). This can also be achieved by directly 
comparing the measured curve R(p), with a set of simulations, and finding the simulation 
that minimize the square of the differences between the measured and simulated curves. 

Note also that the use of the reflectance intensity R(p) and the slope of InR(p) is 
equivalent to the use of two intensities R(p,) and R(pi), measured at two close distances 
Pi and P2, respectively. 
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Method 2. 

The inverse problem can also be performed, considering properties ofR(p) expressed by 
Equ.(1.8). 

Indeed, it can be derived from Equ.(1.9) that die quantity does not depend on the 
absorption coefficient (i*. 

~-^R{ P> /i t ,/i a ,y) (1 . 12 ) 



This property is confirmed in Fig.9 , where the quantity ^— ^ (a/4 > evaluated 

mm, is plotted as a function of y^* for y=l, 1.9 and 2.5 and reduced albedo a =l, 
0.95 and 0.9 (note that any other choice of p is possible, as long as p is of the order of 1 
mfp»). 

In Fig.9. , the parameter -^-^[p.p^p^y) clearly depends on ju*' and y In contrast, 
the dependence on a ' is almost negligible. 

Therefore, ^and/// can be derived from the parameter — J^(p,r/^ i p ffl y) , calculated 

dp " e / 

from the experimental reflecrance R(p) (see section 5.2). Simultaneous determination of y 
and ///require values of —^R(p,^r t ,p a ,y) at different distances p (at least two). If one 

of the parameters /or p s ' is already known, the determination of the other one can be 
obtained from the value of -~^R{p, p,, fi a ,y) at a single distance. 
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Convenien^analytical appr oximation of ^ [p^ t . A, . r) can be obtained by fitting 



Monte Carlo results to polynomial functions. For example, in the case of fixed 7^1.9, we 
obtained: 

Mi '-0.9 1 62-S2.89x+l 795x 2 - 1 8 1 55x 3 +65428x 4 (1.13) 



where x = J Mp, 21 P = ^ x expressed in [mm" 2 ] and n s ' is 

2 dp v 

expressed in [mm" 1 ]. 

This results is valid for a probe of refractive index of 1.5 and a sample of refractive index 
of 1 .4, optical fibers of diameter 200 jim, NA=0.37 (source and collection). 

Once ju s ' and ^are calculated from the procedure explained above, /s a is calculated from 
the absolute value of R(p) t which highly depends on For given /u s \ ^and p values, the 
dependence of reflectance on fi a is obtained h ^Monte Carlo^ results. From Equ.(1.8) we 
have: 



where and are functions given by Monte Carlo simulation. Particularly, they can be well 
approximated by polynomial functions. For example, for 7*1.9, probe refractive index of 
1.5, sample refractive index of 1.4, optical fibers of diameter 200 |im, NA*=0J7 (source 
and collection): 

h « -0.002257 - 8.171 y + 268.8 y 2 (1.15) 
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/- 0.0 1 3 i 1 +0.05 1 84 Mj '-0.0 1974 & ' 2 + 

0.003217 /t/ J -0.0001992//, w (1.16) 

where y = [^(p.ft.M^r) -/] [mm' 1 ]. //, is expressed in [am' 1 ]. 

Method 3 

Equ.(l,8) also shows that relative measurements of Me> i.e. variation of^i* from a known 
value Mae, is possible by measuring the variation of the parameter ^R^p, : 

V*(A#.A.o»rW^(A/C.>W) =B(/i ait tr,)-B(p.,fif,) (1.17) 

This relation is illustrated in Fig. 10. m the case of /w*0, p=l mm, ;r=1.0, 1.9 and 2.5, 
and for a' = 1 to 0.83. confirms that the influence of y is weak in the 
quantity J R(p, u M ,jJ B ^y) - v , ir t ,M„Y) ■ For known the function ate'; allows 
for. the determination of a relative absorption change 4u a -M<t- Mao- Fig. 10. illustrates the 
case Mao "0, but any other value of Mao is possible. The interesting point is that B(a 'J does 
not depend on the phase function. 
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